Characterization of anomalous Zeeman 
patterns in complex atomic spectra 



Jean- Christophe PaiiS and Franck Gilleron 
CEA, DAM, DIF, F-91297 Arpajon, France 



Abstract 

f~i , The modeling of complex atomic spectra is a difficult task, due to 

the huge number of levels and lines involved. In the presence of a mag- 
netic field, the computation becomes even more difiicult. The anoma- 
lous Zeeman pattern is a superposition of many absorption or emission 
profiles with different Zeeman relative strengths, shifts, widths, asym- 
metries and sharpnesses. We propose a statistical approach to study 
the effect of a magnetic field on the broadening of spectral lines and 
transition arrays in atomic spectra. In this model, the a and tt profiles 
are described using the moments of the Zeeman components, which 
(-H ■ depend on quantum numbers and Lande factors. A graphical calcula- 

^ ^ tion of these moments, together with a statistical modeling of Zeeman 

profiles as expansions in terms of Hermite polynomials are presented. 
It is shown that the procedure is more efficient, in terms of conver- 
gence and validity range, than the Taylor-series expansion in powers 
of the magnetic field which was suggested in the past. Finally, a sim- 
ple approximate method to estimate the contribution of a magnetic 
field to the width of transition arrays is proposed. It relies on our re- 
cently published recursive technique for the numbering of LS-terms of 
an arbitrary configuration. 

1 Introduction 

In astrophysics, the observation of a splitting of spectral lines in the visible 
and UV ranges for a few white dwarfs [T] confirmed the existence of intense 
magnetic fields (0.1 - 10^ MG) as predicted by Blackett [2j. The influence 
of a magnetic field on an atom modifies its emission or absorption lines. 
Thanks to this property, known as Zeeman effect, the detection of magnetic 
fields is possible at large distances, through the measured radiation. The 
linear and quadratic Zeeman effects O S] explain the separation of spectral 
lines and enable one to determine a value of the magnetic field. In the same 
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way, pulsars and neutron stars having an even more intense magnetic field 
(10^ - 10^ MG) have been discovered through their spectrum in the range 
of radio-frequencies and X-rays. There are numerous astrophysical appli- 
cations, either direct or indirect, and requiring sometimes a sophisticated 
theoretical modeling. The methods differ according to the nature of the 
objects studied (see tabled]), the magnitude and the geometry of the mag- 
netic fields, and to the quality of the observation in terms of sensitivity and 
spectral resolution. Moreover, the variations of the magnetic field of stars 
during their rotation bring some information about their global geometry. 
The "spectro-polarimetric" methods exploit the additional recording of the 
circular polarization with respect to the wavelength. This enables one to 
obtain a detailed map of the field [5] through a separation of its components 
parallel or perpendicular to the line of sight. 

Strong magnetic fields are also encountered, for instance, in magneto- 
inertial fusion [6j . Inserting a magnetic field into inertial-confinement-fusion 
capsules before compressing them presents the advantages to suppress 
the electron thermal-conduction losses and to better control the a-particle 
energy deposition. The magnetic fields generated inside a Hohlraum can 
reach a few MG. 



Magnetic field B (MG) 


Astrophysical object 


10^ - 10« 
10-1 - 10^ 
10-"^ - 10-2 

- 10-6 
10-13 - 10-11 


Neutron star or pulsar 
White dwarf 
Hot magnetic star 
Planets of the solar system 
Interstellar cloud 



Table 1: Orders of magnitude of magnetic fields encountered in astrophysics 
(1 MG=106 G=100 T). 

In this work, the effect of a magnetic field on the broadening of spec- 
tral lines and transition arrays in complex atomic spectra is investigated. 
A proper description of physical broadening mechanisms |8| requires a si- 
multaneous treatment of Stark and Zeeman effects, which was performed by 
Ferri et al. [9] in the framework of the Frequency Fluctuation Model |10| . 
In the case of an atom (ion) having several open sub-shells, the number of 
electric dipolar lines can be immense and the anomalous Zeeman pattern 
is a superposition of many profiles. When dealing with a huge number of 
simultaneously recorded profiles, it becomes necessary to characterize the 
line shape in terms of a limited number of parameters, and therefore to de- 
termine constraints on modelings. A statistical analysis can be performed 
using the moments of the profile. The n*'^-order centered moment /i„^c[^] of 
a distribution A{E) is defined by 
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JZ,iE-f,,rAiE) dE 
I-ooME) dE 



(1) 



where 



JZ E A{E) dE 



(2) 



IZME) dE 



is the center of gravity of A{E). Each absorption or emission profile con- 
stituting the anomalous Zeeman pattern has its own strength, shift (first- 
order moment), width (second-order moment), asymmetry (third-order mo- 
ment) and sharpness (fourth-order moment). We discuss different ways of 
calculating these moments (whatever the order) in terms of the quantum 
numbers and Lande factors of the levels involved in the line and present a 
statistical modeling of the Zeeman profile. It relies on the use of a A-type 
Gram-Charlier expansion series for each of the components AM=0, +1 and 
— 1. Finally, leaning on our recently published recursive approach for the 
numbering of LS-terms of an arbitrary configuration pT] , we propose a sim- 
ple approximation to estimate the contribution of a magnetic field to the 
emission and absorption coefficients. 

The paper is organized as follows. In section [51 the intensity distribu- 
tion of an electric-dipolar (El) line is introduced, together with its strength- 
weighted moments. In section [3l a graphical representation of the angular- 
momentum sum rules involved in the calculations of the moments is de- 
scribed. It reveals the way the Racah algebra proceeds and is simple to 
compute: the re*''-order moment reduces to a regular polygon with (n -|- 2) 
sides. In section [H the statistical modeling of a line perturbed by a mag- 
netic field is discussed, using particular distributions involving the reduced 
centered moments of the Zeeman vr and a± components. It is proven that 
the Gram-Charlier development is more efficient than the usual Taylor-series 
expansion. In section O an efficient approach to take into account the ef- 
fect of a magnetic field on a transition array is proposed. In section [U] it is 
shown that the techniques presented in this paper still apply when hyperfine 
interaction is included and section [7] is the conclusion. 

2 Intensities and characteristics of Zeeman com- 
ponents 

The Zeeman Hamiltonian reads: 



where B is the magnitude of the magnetic field along the z-axis B = 
B Uz, fJ-B the Bohr magneton, gg = 2.0023192 is the anomalous gyromagnetic 



Hz = fisB {Lz + gsSz) 



(3) 
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ratio for the electron spin, and and Sz respectively the projections of 
total orbital and spin angular momenta of the system. For sufficiently weak 
values of the field the off-diagonal matrix elements of Hz that connect 
basis states of different values of J (modulus of the total angular-momentum 
of the system J = L + S) will be negligible compared to the contributions 
of the Coulomb and spin-orbit interactions to the energy. It becomes then 
reasonable to neglect the mixing of basis states of different values of J. The 
energy matrix breaks down into blocks according to the value of J (as in 
the field-free case) and the contribution of the magnetic field to the energy 
can be calculated as a simple perturbation. The following expression for the 
diagonal matrix element of Hz for the state I7JM) 

(7JMIL, + gsSz\-fJM) = g^j{jJM\Jz\jJM) = g^jM, (4) 

where Jz = Lz+Sz, defines the Lande factor g-yj of level 7 J [12]. One can 
roughly consider that Zeeman approach is no longer valid when the magnetic 
field is of the same order of magnitude as the spin-orbit contribution (see 
table E]): 

Be = {Z*e'^/hc)^me^/{nBh'^). (5) 
In that case, a Paschen-Back [T3] treatment is necessary. 



Element 


Bc{MG) 


H (Z=l) 


0.0078 


Al (Z=13) 


1.30 


Ni (Z=28) 


6.10 


Nb (Z=41) 


13.10 


Sm (Z=62) 


30.00 


Po (Z=84) 


55.00 


Np (Z=93) 


67.50 



Table 2: Critical value of the magnetic field for which the spin-orbit inter- 
action becomes of the same order of magnitude as the magnetic interaction. 
This gives an estimate of the critical field for which Paschen-Back treatment 
is more appropriate than Zeeman description. 

In the presence of a magnetic field, the total intensity of transition 
7JM — )• j'J'M' at the energy E reads: 

^i^) = ^ SyjM,yj'M' 

-yJM-^yj'M' 

'X^'yJM,-y'J'M'iE — EyjM,'y'J'M'), (6) 
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where 



E^jM,-y'j'M' = E^j^^iji + ^bB (gyj'M' — g^jM) (7) 

and SyjM,'y'j'M' are respectively the energy and the strength of a tran- 
sition 7JM — )• 7' J'M'. E^j^yj' represents the energy of the hne 7 J — >• 7' J': 

E^j^yj, = Eyj, - E^j = {iJ'\H\iJ') - {jJ\H\jJ), (8) 

where H is the Hamiltonian of the system. The normaUzed profile 
^'yjM,'y'j'M'{E) takes into account the broadening of the line due to ra- 
diative decay, Doppler effect, ionic Stark effect, electron collisions, etc. 

Assuming that the optical media is passive {e.g. there is no Faraday 
rotation), the intensity, detected with an angle of observation 9, is given by 

MM- 

Ie{E) = mE)cos\e) + I^iE)sm\e), (9) 
where the longitudinal intensity is 

^{E) = ^{I+^{E)+I^^{E)) (10) 
and the transverse intensity 

I^iE) = ^ (/+i(i?) +/_i(i?) +2/o(S)) . (11) 
Ie{E) can be written in the form 



+5^/o(£)^ (12) 

Each line 7J — )• 7' J' can be represented as a sum of three helical com- 
ponents associated to the selection rules M' =M -\- where the polarization 
q is equal to for vr components and to ±1 for a± components. The in- 
tensity of the q component of the El line 7J — )• 7' J' reads, assuming that 
all quantum states are populated in the statistical-weight approximation 
(high-temperature limit): 



h(E) — ^ SM,M',q 
M,M' 

x^jjM,YJ'M'{E — E^jM,YJ'M'), (13) 



where 
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and 



The quantity 



SM,M',q = CM,M',q X Sjjyji (14) 



CM,M',q - 3 ( }^ . (15) 



S'YJ,'y'J' — ^ S^JM,YJ'M' (16) 

■yJM-^yj'M' 

represents the strength of the Hne 7J — )• 7' J' and is proportional to 
|(7J|Z|7' J')p, where Z is the z component of the dipole transition operator. 
Since 

5^ ( -M -q M'J "3' ^^'^^ 

M,M' 1 / 



each component has the same strength. The number of transitions in 
each component is equal to 2x min(J, J')+l. The distribution Iq{E) can be 
characterized by the moments centered in E^j^yj': 



M,M' 
M,M' 

1=0 ^ ^ 

where 

■A^a = E ^MM',, {gj'J'M' - g^jM - M^^y (19) 

M,M' 

is the n*'*-order centered moment of the distribution. It is useful to 
introduce the reduced centered moments defined by 

«n - 2^ CM,M',q [ | , (20) 



M,M' 
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where Ai f is the center-of-gravity of the strength-weighted component 

energies (relative to E^j^yji and in units of hbB) and VV^ = \J -^2^ is 
the standard deviation (in units of ^bB). Centered moments of cr_ and (T+ 
components are related by A^l^^c' = (— 1)"'-A^n^c'- The use of a^n instead 
of M^i" allows one to avoid numerical problems due to the occurence of 
large numbers. The first values are Qq^' = 1, aj^' = and o}-^ = 1. The 
distribution Iq{E) is therefore fully characterized by the values of A^J^', V^''' 
and of the high-order moments with n > 2. It is reasonable to consider 
that the first four moments are sufficient to capture the global shape of 
the distribution Iq{E) (see for instance Ref. [16], p. 88-89). The third- and 
fourth-order reduced centered moments Og^' and af'^ are named skewness 
and kurtosis. They quantify respectively the asymmetry and sharpness of 
the distribution. The kurtosis is usually compared to the value a^^^ = 3 for 
a Gaussian. 



3 Moments of the Zeeman components tt, cr+ and 
cr of a line 7 J — > 7' J' 



3.1 Racah algebra and graphical representation 

The moments can be easily derived using Racah algebra and graphical 
techniques [13 HSl HH [20l [21]. We define the notations [a,6,c, •••] = 
(2a + l){2b + l)(2c + 1) • • • , and use the convention of Biedenharn et al: 
X = x{x + 1) [2^. Since M can be expressed as 



M = i-iy-\/[j]j{ ; ^ ), (21) 



the first-order moment can be obtained from the relations ()8ip . ()83p and 
given in appendix A [T7j. One has 



and 



-M -q M' / 12 

M,M' ^ 1 / 



E(_m M')'«' = ¥''-'^'^' (23) 



M,M' 

which gives finally [23 



M^f = ^ [2{g^j+gYj') 



4 

+(57J - 9yj'){J - J'KJ + J' + 1)] • (24) 
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Figure 1: Graphical representation of a three-jm coefficient. 




Figure 2: Graphical representation of a six-j coefficient. 



The variance is obtained using the sum rule (|82p |17] together with the 
expressions ([83]) to dlH]) [2111251 [26]. More generally, the n*''-order moment 
involves the following sum rule: 

M,M' ^ -i / \ / 

where n is an integer. Figures [1] and [2] give the graphical simplified 
representations of a three-jm and a six-j symbol respectively. Each line 
represents an angular momentum |17[ \TE[ \T9\ \20\ [2T] . The names of the an- 
gular momenta (or of their projections in the case of three-jm coefficients) 
and the phase factors are omitted. Figures [SlHl [5] and [6] display the graphical 
representations of the calculations of the first four moments A4^^^ and 
A^3^' and A4^^^ respectively. One can also see on Fig. [3] how the three-jm 
symbols merge into a single closed diagram. These schemes are a representa- 
tion of summation rules and reduction formulas. Although some computer 
programs exist (see for instance [271 HSl HH [301 [Ml [32] ) , which are devoted 
to the reduction of graphs, it is easy to understand that the calculation be- 
comes more and more cumbersome as the order of the moment increases. 
The n*^-order moment reduces graphically to a polygone with {n+2) sides. 
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Figure 3: Simplified graphical representation of the sum rule ([8T]l of ap- 
pendix A involved in the calculation of the first-order moment A4^^^ of a 
Zeeman component. The first equality corresponds to the merging of the 
three three-jm symbols, and the second one the splitting into a three-jm 
and a six-j symbol. 




Figure 4: Simplified graphical representation of the sum rule ()82p of ap- 
pendix A involved in the calculation of the second-order moment A^2^' of a 
Zeeman component. 
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Figure 5: Simplified graphical representation of the sum rule involved in the 
calculation of the third-order moment M^^^ of a Zeeman component. 




Figure 6: Simplified graphical representation of the sum rule involved in the 
calculation of the fourth-order moment Ai^^^ of a Zeeman component. 
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3.2 Expression in terms of Bernoulli polynomials 



Mathys and Stenflo \33\ [3l] have obtained more compact formulae for the 
moments in terms of Bernoulh polynomials (see appendix B). Values of 
and 04 for the three selection rules A J = 0, —1, +1 are displayed in tables 
[3] and [H One finds that the variance of the tt component is always larger 
than the variance of the 0"+ and cj_ components, indeed: 

yW _ y[-±] = (g7J-gyjO' (gj _ 1) if / = J (26) 

and 

vW_V[a±]^ (glJ^g7Zl!(j + j) if j' = j+i, (27) 
20 

where V'^l = M^^^ = M^^^ - Therefore, in ah cases, - 

yW±] > 0. We can see on Fig. [7] that the variance of the vr component 
for a given value of J is larger for AJ = than for AJ = ±1 lines, and 
that the difference increases with J. Things are slightly different for the a 
components (see Fig. [7]): the variance for A J = overcomes the one from 
AJ = ±1 only for J > 3. Moreover, the difference between both variances 
at fixed J is smaller than for the vr component. Figure [8] shows that the 
skewness of the 0"+ component is a decreasing function of J for A J = ±1 
line (the skewness is zero for AJ = since the splitting is symmetric in 
that case). On the contrary to the variance, the kurtosis 04 (see Fig. E]) 
is systematically higher for AJ = ±1 than for AJ = 0, and the difference 
is almost constant and equal to 1. It is interesting to plot 04 versus Q3 for 
the 0"+ component; it reveals that the dependence is quite linear, and that 
the values are very concentrated around 0.875 for the kurtosis and slightly 
above 3 for the skewness (see Fig. [TO]l . As can be shown on Fig. \TT\ for a 
given value of J the reduced centered moments a„ increase with the order 
n, and, for a given value of n, they increase as well with J, and get closer 
and closer when J increases. 



0-+ 


J' = J 


J' = J + 1 


J' = J - 1 


as 
a4 




5 / 12j-17\ 
7 1, 4J-3 J 


2V5 J+1 
3V3 ^/.J(J+2) 
5 /l3J(J+2)-4\ 
21 1^ J(J+2) ) 


2V5 J 
3V3 VJ2-1 

5 f 13J2-17^ 
21 1-J^ ) 



Table 3: Values of as and a4 of the a+ component of El lines. 

The numerical values 04, ag and ag of the vr component for several 
lines are listed in table O Tables [6] and [7] contain the odd reduced centered 
moments of the a component for the same lines. 
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TT 


J' = J 


J' = J+l 


J' = J - 1 


as 
a4 




25 /3{(J+2)j2-l}J+l\ 




5 / 3J(J+2)-2\ 




5 / 3j2-5\ 


7 1, {i-3j}2 ; 


7 V j(j+2) ; 


7 i j 



Table 4: Values of as and of the vr component of El lines. 




Figure 7: (Color online) Variance of the vr and a components with respect 
to J. 

4 Zeeman profile in low magnetic fields 

In the following, we consider the case where 

I{E)='^{I+^{E) + I^^{E)+h{E)) (28) 

which, according to Eq. ([9]), corresponds to an observation angle 6 with 
z axis such that cos^(0) = ^. 

4.1 Taylor-series expansion 

In the following, we make the assumption that ^-yjM,YJ'M' is a universal 
function centered in E^jM,'y'j'M' ■ The quantity Iq{E) (fT3|l can be ex- 
pressed [331 [Mj ^-s a Taylor series around the line energy E^j^^ijr. 
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J 

Figure 8: Skewness 03 of the 0"+ component with respect to J for A J = ±1. 



Line 


V 


04 


"6 


as 




0.60 


1.667 


2.778 


4.629 




0.03 


2.048 


5.190 


14.407 


"^-03/2 ^-^5/2 


1.05 


1.871 


3.944 


8.436 


5P2^ 


1.60 


1.964 


4.576 


11.230 



Table 5: Parameters of the vr component for several lines of the transition 
array Fe VII 3d^ — )• ?>dAp. Even reduced centered moments. 



fc=i 

x^*(^-^7W)|- (29) 
Assuming a Gaussian physical broadening of the lines: 
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3 



3.2 
2.8 
2.4 
2 
1.6 
1.2 



1 ' 1 ' 1 ' 1 


' 1 










✓ 

✓ 


— 7i: J'=J+1 or J-1 




f 


■-■ Ji: J'=J 




1 

1 


-- o: J'=J+1 or J-1 




1 

I 


- ■ o: J'=J 




1 






/ /' 






/ / 






1 f' 

/ / 
/ / 

/ / 






" 1 1 , 1 , 1 , 1 , 1 " 
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Figure 9: (Color online) Kurtosis 04 of the vr and a components with respect 
to J for AJ=0 and AJ = ±1. 

where v represents the variance of the physical broadening mechanisms 
other than Zeeman effect (Doppler, Stark,...), we have (Rodrigues' formula): 

(-1)" 

xHe„ (^^^^) , (31) 

where He^ is the Hermite polynomial of order related to the usual 
Hermite polynomial by 

= (^) • (32) 

He/j obeys the recursion relation 

Hefc+i(x) = X Hefc(x) - A; Hefc_i(x) (33) 
with Heo(x)=l and Hei(x) = x. The resulting expression of ^q{E) reads 
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1 1.05 
Skewness a. 



Figure 10: Kurtosis 04 versus skewness 03 of the 0"+ component for A J=+l. 



At the second order 



k=l 



(34) 



IqiE) — S^j^YjiX'^{E — E^j^Yji) 



X < 



1 + 



[g] (E - E^j^Yj> 



-{y^BBf (M'!i + (M' 



V- (E-E, 



■yJ,YJ' 



2^2 




(35) 



Throughout the paper, the calculations denoted "exact" are performed 
with the Flexible Atomic Code (FAC) code fSSj- Figure [12] shows that, for 
i?=1.25 MG and v=5 10~^ eV^, the TS expansion converges to the exact 
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J I \ I I I \ I L 

3 4 5 6 7 

n 



Figure 11: (Color online) Reduced centered moments a„ of the C7_|- compo- 
nent versus n for different values of J in the case A J = ±1. 

profile with a very good accuracy. Such an approach still works fairly well 
even when the profile starts to exhibit oscillations due to the important 
separation of the vr, a+ and cj_ components (see Fig. [T3|l for B=1.5 MG 
(corresponding to fj,BB/y/v ~ 1.23). In the latter case however, the conver- 
gence is quite slow: a satisfactory agreement is still not achieved at the order 
n=16. The Taylor-series method is valid for fisB < ^/v, but breaks down 
if fisB becomes much larger than ^/v. Note that expression (|35p can be 
exploited for a rough determination of the magnitude of the magnetic field 
B, provided that variance v of the other broadening mechanisms is known 
(see appendix C). It is interesting to mention, as can be seen on Fig. [T^ 
that the modeling of each component separately is not satisfactory at all, 
since in the present case, the separate TS expansions exhibit some oscilla- 
tions and can even become negative, for (T+ and (T_ components. However, 
such variations do not affect the resulting total function (sum of the three 
components). 

4.2 A-type Gram-Charlier expansion series 

An alternative to the Taylor-series method consists in using a statistical 
distinction based on the Gram-Charlier development. Once the centered 
moments fj,n,c of a discrete distribution A{E) are known, such a distribu- 
tion can be modeled using an analytical function which preserves an arbi- 
trary number of these moments. It is possible to build a function using 
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Figure 12: (Color online) Modeling of a line J = 3 — )• J' = 4 of transition 
array Fe VII 3(P 3d4p with Taylor -series expansion of different orders, 
compared to the exact calculation for B=1.25 MG and v=5 10~^ eV^. The 
development for n=16 and the exact calculation are superimposed. 

the properties of orthogonal polynomials and their associated basis func- 
tions [m [Ml EZl [38] . The A-type Gram-Charlier (GC) expansion series is a 
combination of products of Hermite polynomials by a Gaussian function: 

expf-^) / n \ 
GCniE) = ^ / (l + Y^Ck Hefc(y) , (36) 

with 

where y = {E — Aii)/^/^2,c) is the number of moments, m.t{k/2) is 
the integer part of fc/2 and the Hermite polynomial He^ is defined in the 
preceding subsection 14.11 The GC series uses the reduced centered moments 
Q-nlM of which are defined by: 
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30 



3 



o 20 



•■^ 10 



— Exact, B=1.5 MG 
TS, n=8 
TS,n=16 




53.18 53.2 53.22 

Energy (eV) 



53.24 



53.26 



Figure 13: (Color online) Modeling of a line J = 3 — )■ J' = 4 of transition 
array Fe VII 3d? — )■ SdAp with Taylor-series expansion of different orders, 
compared to the exact calculation for B=1.5 MG and v=5 10"^ eV^. 



The fourth-order GC series reads: 



+^^(3-6n^ + n^)}. (39) 

The truncated series GCn{E) may be viewed as a Gaussian function 
multiplied by a polynomial which accounts for the effects of departure from 
normality. Therefore it may be a slowly converging series when A{E) differs 
strongly from the Gaussian distribution. It is also known to suffer from nu- 
merical instability since Eq. (|37|) involves a sum of large terms of alternating 
sign. Still assuming a Gaussian physical broadening (see Eq. ([30]) 1 of the 
lines, the moments of the convolution read: 
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I I I I I I I I I 

53.18 53.2 53.22 53.24 53.26 53.28 
Energy (eV) 



Figure 14: (Color online) Modeling of a line J = 3 — )• J' = 4 of transition 
array Fe VII 3(P — )• 3d4p with Taylor-series expansion of eighth order, for 
each component separately for B=1.5 MG and v=5 10~^ eV^. 

/I + (-l)n-k\ 

x[ \' yk,M], (40) 

where x i— )■ r(x) is the usual Gamma function. 

4.2.1 Global Gram-Charlier expansion series for the total inten- 
sity 

In that case, A = I = J2ql-i 
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Line 


"3 


as 






0.994 


5.521 


27.913 




0.889 


6.067 


41.822 


^1)3/2 ^ '^5/2 


0.939 


5.856 


35.177 


5P2^ 


0.913 


5.977 


38.670 



Table 6: Parameters of the a component for several lines of the transition 
array Fe VII 2>d? — )• 3(i4p. Odd reduced centered moments. 



Line 


V 


04 


"6 


as 




0.4500 


2.778 


12.654 


62.592 




2.2500 


3.032 


16.426 


114.19 


4Z)3/2 ^ 4D5/2 


0.7875 


2.914 


14.637 


87.850 


5P2^ 


1.2000 


2.976 


15.575 


101.273 



Table 7: Parameters of the a component for several lines of the transition 
array Fe VII 3d^ — )• 3(i4p. Even reduced centered moments. 

i=0 j=0 ^ J \-> J 

Figure [15] shows that, for a line J = 3 — )• J' = 4 of transition array Fe VII 
3d? — )• 3d4p, the fourth-order A- type Gram-Charlier distribution GC/i{E) of 
Eq. (j39|) provides a satisfactory depiction of the profile. However, when 
the order increases, the departure from the exact calculation becomes larger 
and larger. This is due to the fact that the reduced centered moments an\I] 
(see Eq. (f38l) ) do not depend on B. Therefore, such an approach can be 
applied only if the global shape I{E) is close to a Gaussian, i.e. does not 
have a non-monotonic character. This implies that the method is valid only 
if hbB < ^/v, so that the vr and a± components are not too separated. This 
approach provides a good depiction of the profile if fisB < ^/v. 

4.2.2 A-type Gram-Charlier expansion series for each compo- 
nent 

In that case, A = Ig and 
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— Exact, B=1.5MG 
GC global, n=4 




Energy (eV) 



Figure 15: (Color online) Modeling of a line J = 3 — ?• J' = 4 of transition 
array Fe VII 3d? — )• SdAp with A-type Gram-Charlier expansion series of dif- 
ferent orders, compared to the exact calculation. B=1.5 MG and v=5 10~^ 
eV^. The sum of three fourth-order GC functions and the exact calculation 
are almost superimposed. 



/^n,c[/,]=M',U(^Bi?r. (42) 

This approach has a wider validity range than the previous one (see 
Fig. [16] the case of a magnetic field equal to B=2.5 MG). When the ratio 
fiBB/^/v becomes larger than one, the summation of three A-type Gram- 
Charlier expansion series brings more flexibility. One can notice on the wings 
that A-type Gram-Charlier expansion series yield negative values in certain 
circumstances. However, it provides a good global depiction of the profile. 
Figure [T7] displays the modeling of each component separately. The 0"+ and 
cr_ profiles do not show the oscillations observed with the TS expansion (see 
Fig. dH). 
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15 



53 



— Exact, B=2.5 MG 
-- GC global, n=4 

— GC global, n=8 

— GC +GC +GC , n=4 



53.1 




53.2 
Energy (eV) 



53.25 



Figure 16: (Color online) Modeling of a line J = 3 — t- J' = 4 of transition 
array Fe VII 3d? — t- SdAp with A-type Gram-Charlier expansion series of 
different orders, compared to the calculation. B=2.5 MG and v=5 10~^ 
eV^ . The sum of three fourth-order GC functions and the exact calculation 
are almost superimposed. 

5 Global accounting for Zeeman effect on a tran- 
sition array 

5.1 Statistical description 

The absorption and emission spectra consist of a huge number of electric- 
dipolar (El) lines. A transition array [39] represents all the El lines between 
two configurations and is characterized by a line-strength-weighted distri- 
bution of photon energy E: 



■yJ^-y'J',M,M',q 

^ S^j^-y'j' 'if{E - E^j^yj/). (43) 

The sum runs over the upper and lower levels of each line belonging to 
the transition array. 
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53.1 53.15 53.2 53.25 

Energy (eV) 



Figure 17: (Color online) Modeling of a line J = 3 — )• J' = 4 of transition 
array Fe VII 3(P — )• SdAp with A-type Gram-Charlier expansion series of 
fourth order for each component. B=2.5 MG and ^=5 10~^ eV^. 

In the UTA (Unresolved Transition Arrays) approach [3D], the discrete 
set of lines (as 6 functions) is replaced by a continuous function (usually 
Gaussian) which preserves its first- and second-order moments. The mo- 
ments of this distribution are evaluated as 

It is possible to derive analytical formulae for the moments /i^ using 
Racah's quantum-mechanical algebra and second-quantization techniques 
of Judd [41J- Such expressions, which depend only on radial integrals, have 
been published by Bauche-Arnoult et al. [401 [32l [33l 133] for the moments 
fJ'n,c (centered moments with respect to ^i) with n < 3 of several kinds 
of transition arrays (relativistic or not). Karazija et al. have proposed 
an algorithm in order to calculate the moments of a transition array using 
diagrammatic techniques [271 EH [29] . 

The contribution of Zeeman effect to the fc^'^-order moment of a transition 
array for a polarization q reads [35] 
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I^k = ^ 'P-yjM^yj'M' 

7JM,7'J'M' 

X {E^j,y,j, + fiBB {gyj'M' - 57JM)}' , (45) 

VyjM^'y'j' M' being the probability of a transition from {J,Mj) to {J' ,Mj — 
q) (component). Using the binomial development, one obtains: 

^^f = l^k + Y^^I^bBY M Hk-pJ'p, (46) 
p=l ^-^ ^ 

where 

l^k= ^ 'P-yJM^yj'M'iEyj/ - E^j)'' (47) 
7JM,7'J'M' 

and 

•^P= {gyj'M' - g^jM)PV^jM^yj'M', (48) 

7JAf,7'J'M' 

which can be evaluated using the techniques mentioned in the preceding 
sections [3] and HI 

The complexity of such a calculation encouraged us to develop an al- 
ternative approximate method. Suppose one wants to include the effet of a 
magnetic field in a numerical code devoted to the computation of opacity 
or emissivity, without performing the diagonalization of the Zeeman Hamil- 
tonian. The numerical code can be either based on a detailed (see sections 
[21 [3] and HD or a statistical description (relying on the UTA formalism as 
mentioned above). The main contribution comes from the splitting of the 
line into three components. Indeed, if one considers 3 components with zero 
width positioned at E^j^yj' — fisB, E^j^y ji and E^j^y ji + hbB (each having 
the same strength S^j^yj'), the variance is equal to: 



S-yJ,YJ' V 3 

3 



"I 5 [E^j^yji - E^j^yji) 



H ^'''^ {E^j^yji + ^bB - E^j^yji)"^^ 



(49) 



which is equal to 2/3 {nBBf ^ 3.35 lO'^ [5(MG)]2. 
The broadening of each q component separately due to the magnetic field 
(which is larger for a tt than for a a component as a consequence of Eqs. 
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(|26p and (|27p ) is always much smaller than 2/3 {^bB)"^ (by at least one 
order of magnitude). Thus, the contribution of a magnetic field to an UTA 
can be taken into account roughly by adding a contribution 2/3 {fisB)'^ to 
the statistical variance. In case of a detailed transition array, the Zeeman 
broadening of a line can be represented by a fourth-order A-type Gram- 
Charlier expansion series (Eq. ()39p ). i.e.: 



^z{E-E^j^yj,) = Yl 



exp[-4- 



q=-l fXBB^2TrM^^l 



where 



^iBB^Ml[ 

The coefficient ge of the line 7 J — )• 7' J' is given by 



(50) 



^ ^ - E^j^yjf -q ge fisB ^^^^ 



9e = ^ {2(57 J +9YJ') 

+ {9^J - 9^'J'){J - J'){J + J' + 1)} , (52) 

where g^j and g^iji are the Lande factors of levels 7J and 7' J' respec- 
tively [MlllZlllH]. 

5.2 Approximation of the coefficient 

If the values of g^j and g^iji are unknown, we suggest to replace ge by its 
average value in LS coupling ge- Knowing the distribution of spectroscopic 
terms Q{S,L) [iQl [11], it is possible to get a quick estimate of ge- Indeed, 
the equality 

XsLJ = Y,QiS^^)^SLJ, (53) 
7(5L)J S,L 

where Xslj is any quantity depending on 5, L and J, enables one to 
deal with the coupling of angular momenta L and S avoiding the use of 
coefficients of fractional parentage. One has 
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9e=Yl Q^^^^'^9e{S,L,J,L',j')e{L,L',J,j'), (54) 

S,L,J L',J' 

where e{L, L' , J, J') stands for the selection rules: L' = L,L — 1 or L + 1 
avoiding L' = L = and J' = J, J — 1 or J + 1 avoiding J' = J = 0. One 
has 

ge{S,L,J,L',J') = ^ {2{gsLj + gsL'j') 

+ {9SLJ - gSL'J') 

x(J-J')(J + / + !)}, (55) 
where the Lande factors are estimated in LS coupling: 



9SLJ 



1 + 



i9s 



9s + I 



1){J + S-L) 

1){S-L) 



+ 



2J 

(.9s - 



2J 



(56) 



with the convention of Biedenharn et al. [22], x = x{x + 1). The quan- 
tity gs represents the anomalous gyromagnetic ratio defined in section [TJ 
Assuming g'^ ~ 2, one has 



9SLJ 



3 

2 + 



2J 



(57) 



Table [8] contains values of the Lande factor calculated in LS coupling 
using Eq. ([57]) as well as geiS, L, J, L' , J') factor for different lines. 



Line 


9SLJ 


9SL'J' 


ge{S,L,J,L',J') 







1 


1.5 


^Dl^ 7^2 


3 


2 


1.5 


4^3/2 ^ ^D,/2 


1.2 


1.371 


1.417 


5P2^ 


1.833 


1.667 


1.5 



Table 8: Lande factors for several lines evaluated from formula (|57p . 

The problem of listing the terms arising in a complex configuration can 
be solved from elementary group theory [50\ \5T\ \52\ [53l . The number 
Q{S, L) of LS terms of a configuration £^ can be obtained from the relation 

5+1 L+1 

Q{S,L)=J2 2] (-l)''-''^+''-'^"^'iv(M5,Mi), (58) 

Ms=SA-Il=L 
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where Pn{o,, b), number of states with a given Ms = a and Mi = b, can 
be obtained using recursive formulas |llj : 



^ G N 

P^iMs,ML) = -Y.Y.(-^)'''' 



1=1 k=l 



xP^.fc ( Ms - k^-^; ML-kmi+^ ' 



2 ' 2 , 

(59) 

where G = EiIiG^,G, being the degeneracy of orbital i. For the non- 
relativistic configuration l-^ : 

n...^^^-"^ ; ^<k<U^2 (60) 

and for the relativistic configuration : 

mk = k-j + l ; 1 < A; < 2j + 1. (61) 
The recurrence (1591) is initialized with 



Po{Ms, Ml) = 5{Ms)5{Ml). (62) 
For a configuration £^^£2^1^'^ • • • , P(Ms, Mi) is determined through the 



relation 



Pn,,n,AMs, Ml) = (Pn, ^Pn,^--- )iMs, Mi), (63) 
where the distributions are convolved two at a time, which means 



{Pn,(S)Pn,){Ms,Mi) = 
+00 +00 

^ ^ Pn.{M's,M'^) 

M'g=-oo Af£=-oo 

X Pn^ {Ms -M's,Mi-M'i). (64) 

Thus, in order to take into account approximately the impact of the 
magnetic field when the number of lines is large, we suggest to convolve the 
transition array in the absence of a magnetic field with the distribution of 
Eq. ([50]). Figures UHl US] and [20] (example taken from Mc Lean [551 156] ) 
show that the results are quite close to the exact calculation. The main 
approximation here comes from the fact that ge is replaced by its average 
value in LS coupling, which is justified in case of very strong magnetic fields 
(Paschen-Back effect) . Table M displays the energies and Lande factors of 
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Energy (eV) 

Figure 18: (Color online) Effect of a 15 MG magnetic field on transition 
array Fe VII 3d? — t- Sd'ip with a convolution width of 0.017 eV. The curves 
corresponding to the exact (-6=15 MG) and statistical calculation are almost 
superimposed. 



the levels of configurations ls2s and ls2p in intermediate coupling and table 
[To] indicates the oscillator strength multiplied by the degeneracy 577/77,7' j' 
of the six lines. 

If needed, the evaluation of ge can be refined. For instance, it is possible 
to calculate an average value of g.yj depending only on J. This can be 
achieved using the sum rule [57] 

^9-(J = J2 9LSJ^{lJWLSjf = ^ QLSJ, (65) 

7 aLS 7J aLS 

which states that the sum of the Lande factors for any given J is inde- 
pendent of the coupling conditions. Such a property stems from the fact 
that the trace of a matrix is invariant under an orthogonal transformation. 
One can thus define an average Lande factor associated to a given value of 
J: 

T.^9^J ^ EaLs9LSJ . . 

Q{J) Q{J) ' ^ ^ 



where Q{J) is the number of levels having angular momentum J 
which can be evaluated recursively [TT], in a similar manner to Q{S, L) (see 
Eq. ([58])). 



28 



0.4 



a 

3 



,0.3 



S 0.2 



a 

H 



00 



0.1 



Exact, B=0 
Exact, B=15 MG 
Statistical 




52.5 



53 53.5 
Energy (eV) 



54 



Figure 19: (Color online) Detail of Fig. [18 



6 Hyperfine structure 

The same methodology can be applied in order to determine analytically the 
moments of the hyperfine components of a line. The hyperfine operator in 
the subspace corresponding to the relevant nucleus and atomic level reads: 



= A J (I. J), 



(67) 



where Aj is the magnetic hyperfine-structure constant of the level 7J. 
The n*'^-order moment of the hyperfine components is provided by the ex- 
pression 



Mn 



[{jJIFMlHmhJIFM) 



Sppi M M' 

ff MM p^p,^,^^]^, 

-{iJ'lF'M'\H„,WJ'lF'M')Y 
x{'yJIFM\Z^^^\-fJIFMf, 



(68) 



where Zg is the g-component of the dipole operator Z^^\ The J-file 
sum rule |58) enables one to simplify the expression of the strength: 



29 




Figure 20: (Color online) Effect of a 1 MG magnetic field on C V ls2s 
ls2p triplet transition with a convolution width of 0.005 eV. 



Sff'mm' = {lJIFM\Z^^^\-fJIFM)' 

F,F',M,M' 

F,F',M,M' \ ^ / 

F,F' 

= ^E[^] = ^[^"^]' (69) 

and therefore 

L ' ' F F' 

x{F\\zW\\F'f, (70) 
where Xfij = (—1)'^ {F — I — J). Equation ([TOjl can be written 
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Level number 


J 


Energy (eV) 


Lande (IC) 


Configuration 


1 





303.99067 


1.500000000 


ls2s 


2 


1 


297.88824 


2.002320051 


ls2s 


3 





303.59204 


1.500000000 


ls2p 


4 


1 


303.59057 


1.501152782 


ls2p 


5 


1 


307.61242 


1.000007244 


ls2p 


6 


2 


303.60607 


1.501160026 


ls2p 



Table 9: Energy (relative to the energy of Is^) and Lande factor of the 
different levels of configurations ls2s and ls2p. 



Initial level 


Final level 


g'yjf-yJ,YJ' 


4 


1 


1.56936 10-^ 


1 


5 


9.47301 10~2 


2 


3 


3.99278 10-2 


2 


4 


0.11969 


2 


5 


3.80306 IQ-^ 


2 


6 


0.20016 



Table 10: El lines in the transition array l,s2s — )• ls2p in intermediate 
coupling. 



or 



Mn = -^Y.i^,F']{AjXFjj-Aj,XF,jj,r 

^ ' ^ F,F' 

x((/J)F||2«||(/J)F')^ (71) 



Mr 



^ [F, F'] (AjXfij - Aj,Xp,ij,) 

F,F' 

F 1 F' 



In the case where F or F' is equal to 0, the calculation is very simple 
[S9] . In the general case, using 

X„j = 2(-l)^+^+yiJ[/,J]| [ (73) 
one has to calculate: 
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F,F' ^ ^ 

F' J' I ( F 1 F' 



ki 



1 I J' ] J' I J 



(74) 



which can be done usmg graphical methods [17j. Another approach 
consists in adopting another point of view, leading to the evaluation of 
quantities of the type: 

Sn = ^[F]{F-ar (75) 

F 

where a is a constant (depending on other quantum numbers). Such a 
quantity can be expressed, as for the Zeeman effect, in terms of Bernoulli 
numbers (see appendix B): 



n k , \ / , 



k=0 j=0 

X {(2(Sfe+,+2(/ + J + 1) - Bk+,+2) 

+Bk+j+i{I + J + 1) - Bk+j+i)/{I + J + 1) 
+{2{Bk+j+2{\I - J| + 1) - Bk+j+2) 
+Bk+,+i{\I - J| + 1) - Bk+,+i)/{\I - J| + 1)} . 

(76) 

The splitting of F components in a weak magnetic field |60] is in ev- 
ery way similar to the splitting of J levels. The scale of the splitting is 
determined by the factor gp, which is defined by 

{jJIFM\H,\jJIFM\) = ^sBg^j^-^^^M, (77) 
and connected with the Lande factor by 

F + J-I 

9F = g-yj ^ • (78) 

7 Conclusion 

In this work, a statistical modeling of electric dipolar lines in the presence of 
an intense magnetic field was proposed. The formalism requires the moments 
of the Zeeman components of a line 7J — >• 7' J', which can be obtained ana- 
lytically in terms of the quantum numbers and Lande factors. It was found 
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that the fourth-order A-type Gram-Charher expansion series provides bet- 
ter results than the usual development in powers of the magnetic field often 
used in radiative-transfer models. Using our recently published recursive 
method for the numbering of LS-terms of an arbitrary configuration, a sim- 
ple approach to estimate the contribution of a magnetic field to the width 
(and higher-order moments) of a transition array of El lines was presented. 
We hope that such results will be useful for the interpretation of Z-pinch 
absorption or emission spectra, for the study of laser-induced magnetic fields 
in inertial- fusion studies, for the modeling of magnetized stars as well as for 
any application involving magnetic fields in spectroscopic studies of atomic 
and molecular systems. 
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8 Appendix A: Expressions involving three-jm and 
six-j symbols used in sections [3] and |4] 

[a, 6, c, • • • ] = (2a + 1)(26 + l)(2c + 1) • • • (79) 

x = x{x + l) 
h = J-J'' + 2 



J 1 J' \ f J' 1 J \ f J 1 J 

-M -q M' J \ -M' q M J \ -M M 

1 1 1 

-q q 





J 






1 





(81) 



J 1 J 



^ V -M M 

MM' 



J' 1 J 

-M' q M 



J" 





J J 1 
J" 1 J 



J 


1 


J' 


M - 


-q 


M' 


J 


1 


J 


-M 





M 




J" 


1 







q 


f J' 


J 


1 


1 J" 


1 


J 



(82) 
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1 I M=(_l)«+i^. (83) 

-q q J ^ ' V6 ^ ^ 



1 J" i\ {-iy"{J" -i){J" + 2) 



y v/(2- J")!(J" + 3)! 
J' J 1 1 



1 J 



vm 



(84) 

(85) 



1 1 J 



{-iy+' , (86) 

2V6[J]J 



{2' f 



3 h{Ii - 1) - 8.7 



-M M 



Vl20[J](2J- 1)(2J + 3)J 
(-1)^ 



(87) 



V[J]J' 

9 Appendix B: Bernoulli polynomials and num- 
bers 

The Bernoulli polynomials can be obtained by successive derivation of a 
generating function: 



(89) 

t=o 



One can write 

Bn{x) = J2(l)Bmx--\ (90) 

fc=0 ^ ^ 



where -Bn(O) is the n-order Bernoulli number, which is non-zero only if 
n is even and which can be obtained from the relation: 



k=0 ^ ^ 

The first Bernoulli polynomials are 

Bo{x) = 1, (92) 

Bi{x) = x-^, (93) 
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and 



B2{x) = — X + 



Bsix) = {x- l)(x - -)x, 



BAx) =x'^- 2x^ + x^ - —. 
^ ^ 30 



The Bernoulli polynomials obey the following identity: 

Bp+iin + 1) - Bp+i{0) 



k=l 



P+1 



(94) 
(95) 

(96) 
(97) 



and the Bernoulli numbers have the explicit Laplace's determinantal 
formula |61| : 



Bn{0) 



1 

2! 



1 





1 

1 



(n-1)! 
(n+1)! nJ. 



1 



1 



1 



(98) 



10 Appendix C: diagnostic of the magnetic field 

Using second-order TS expansion p5p and assuming the knowledge of the 
variance v of the other broadening mechanisms, it becomes possible to esti- 
mate the magnitude of the magnetic field from the measurement of the full 
width at half maximum (FWHM) of the line 6 = FWHM/(2/u). 



B 



1 



^iB^ C{e) (1 - 2e-'5V2(i _ ^2)) 



where 



with 



and 



c{e) = A{e) (M'^+])Va^5;+]+d(0)a^w 



A{e) 



l + cos2(0) 



(99) 



(100) 



(101) 



sin2(6') 



(102) 
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This simple formula ()99p can provide an estimation of the magnetic field, 
even if the other broadening mechanisms (Stark, electron collisions, Doppler, 
autoionization) are dominant. However, it is not as efficient as the method 
proposed by Stambulchik et al. [62j, which is applicable in situations where 
the magnetic field has various directions and amplitudes (or if they vary in 
time). 
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